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Abstract 

We consider a hyperbolic-parabolic model of vasculogenesis in the multidimensional case. For this 
system we show the global existence of smooth solutions to the Cauchy problem, using suitable energy 
estimates. Since this model does not enter in the classical framework of dissipative problems, we analyze 
it combining the features of the hyperbolic and the parabolic parts. Moreover we study the asymptotic 
behavior of those solutions showing their decay rates by means of detailed analysis of the Green function 
for the linearized problem. 

1 Introduction 

In this work we present some analytical results on the PDEs model of vasculogenesis proposed by Gamba 
elal. [15]. 

Vasculogenesis is the process of blood vessel formation occurring by the production of endothelial cells 
and is lead by a chemotactic phenomenon, i.e. cells direct their movement according to certain chemicals 
in their environment [4, 18]. At first, it was believed to occur only during embryologic development but 
recently it was realized that vasculogenesis can also occur in the adult organism. Circulating endothelial 
progenitor cells were identified and it was observed that they were able to contribute to neovascularization, 
such as during tumor growth, or to the revascularization process following a trauma, e.g., after cardiac 
ischemia. 

In this paper we proceed in the mathematical study of the following system proposed by Gamba et al. [15] 

d t p + V-{pu) = 0, 

< d t [pu) + V- {pu® u) + V_P(p) = -apu + p,pV<p, (1) 
df(f>- DA(p + ap- b(p. 

Here p is the endothelial cells density, u the cells velocity and <p the concentration of the chemoattractant. 
Moreover, the positive constants D, a, and b are, respectively, the diffusion coefficient, the rate of release 
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and the inverse of the characteristic degradation time of chemoattractant. The other two positive constants, 
a and ji, measure respectively the friction of the cells on the substrate and the strenght of the cell response 
to the chemical signals. 

Classical arguments for hyperbolic-parabolic systems fail in the case of system (1) due to the production 
term in the third equation, i.e.: a > 0. The treatment of these terms is not an easy task so, as a first step in 
the analytical study of this system, we focus on a particular case. Actually, we restrict ourselves to a specific 
class of solutions, i.e.: to small perturbations of non null constant states. Indeed, we are able to show 
the existance of global solutions for this Cauchy problem in the full multidimensional setting, and some 
precise time decay estimates for those solutions. A first mathematical results on system (1) was proved 
by Kowalczyk et al. in [25] , where a viscous term y V 2 v in the second equation is considered to introduce 
an energy mechanism that models the slowing down of cells in the proximity of network structure. They 
performed a detailed linear stability analysis of the model in the two dimensional case, with the aim of 
checking their potential for structure formation starting from initial data which represent a continuum cell 
monolayer. This model is unstable at low cell densities, while pressure stabilizes it at high densities. 
In [7], Di Francesco andDonatelli dealt with diffusive relaxation limits of system (1) toward Keller- Segel type 
systems, either hyperbolic-parabolic or hyperbolic-elliptic. In order to produce a nontrivial class of solu- 
tions to the hyperbolic system which after a proper rescaling relax toward a Keller-Segel type model, they 
provided by means of Friedrich's symmetrization technique and by linearization arguments, an existence 
theorem (local in time) for the approximating system. Moreover they proved the uniform estimates needed 
to justify the assumptions in case of initial densities which are small perturbation of an arbitrary non zero 
constant state. 

More in general, hyperbolic-parabolic systems have been widely studied by Kawa-shima and Shizuta [22, 
23, 34] . Under the smallness assumption on the initial data and the dissipation condition on the linearized 
system, they were able to prove global (in time) existence and asymptotic stability of smooth solutions to 
the initial value problem for a general class of symmetric hyperbolic-parabolic systems. 
However system (1) does not enter in this framework. As a matter of fact, due to the presence of the source 
terms ap, the dissipative condition fails. 

If we linearize the differential part of (1) we get the semilinear hyperbolic-parabolic model, introduced by 
Hillen to describe chemosensitive movements [19], 

d t p + V-v = 0, 

< d t v + j 2 Vp = -c{(p,V(p)v+h((p,V(l))g{p), (2) 
d t (p = A0+/(p,(/)), 

where p is the density of the population with finite speed j,v-pu the flux and <p the chemoattractant con- 
centration. With reference to the one dimensional case, a first result of local and global existence for weak 
solutions, under the assumption of turning rate's boundness, was proved in [21]. Recently Guarguaglini et 
al. in [16] proved more general results of this model under weaker hypotheses, by showing a general re- 
sult for global stability for a zero constant state in the Cauchy problem and for a small constant state in 
the Neumann problem. These results have been obtained using the general theory of linearized operators, 
and an accurate analysis of their nonlinear perturbations. Proceeding along these lines, in [10] the authors 
presented a global existence theorem and the asymptotic behavior for smooth solutions with small initial 
data to the Cauchy problem, for a simplified version of system (2) in the two dimensional case. Moreover 
in [8, 9] it has been considered the multidimensional model (2), and showed the global existence of smooth 
solutions with small initial data to the Cauchy problem and determinated their asymptotic behavior. 
Since in our case it is not possible to apply the technique of the semilinear case, we follow a different ap- 
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proach. The basic idea is to consider the hyperbolic and parabolic equation "separately", and to take ad- 
vantage of their respective properties. Let us explain our approach in more details. 

We look at the hyperbolic part of system (1) without the source term /ipV0, i.e. we consider isentropic Euler 
equations with damping. This system enters in the general framework proposed by Hanouzet and Natalini. 
In [17], they determined sufficient conditions which guarantee the global existence in time of smooth solu- 
tions for small initial data, which are the entropy dissipation and the Shizuta-Kawashima conditions. 
The first one is a condition for systems which are endowed with a strictly convex entropy. Even if the strict 
convexity guarantees that the entropy estimates are equivalent to the L 2 estimate, and the dissipation yields 
the invariance in the same norm, this condition is too weak to prevent the formation of singularities. In- 
deed there exist systems that satisfy this condition such that there is no global solutions for some arbitrarily 
small initial data. The condition (SK) is an adaptation to hyperbolic problems of the Kawashima condition 
for hyperbolic-parabolic ones. In terms of stability it guarantees the necessary coupling between conserved 
and non conserved quantities in order to have dissipation effects, in both the sets of state variables. 
Following the work by Hanouzet and Natalini [17] and Yong [37], our approach is based on energy esti- 
mates for the parabolic and hyperbolic equations. As a matter of fact, even if classical arguments fail in 
the estimate of the source term fipS/cp, we are able to treat it thanks to particular estimates of the parabolic 
equation. 

Once that the global existence for smooth solutions has been obtained for perturbation of small constant 
states, we are able to determine the asymptotic behavior for large times of solutions, by using the decay 
rates of the Green functions. Our strategy consists in using the decomposition of the Green function of dis- 
sipative hyperbolic systems done by Bianchini at al. [2] and its precise decay rates. Indeed in [2] the authors 
proposed a detailed description of the multidimensional Green function for a class of partially dissipative 
systems. They analyzed the behavior of the Green function for the linearized problem, decomposing it into 
two main terms. The first term is the diffusive one, and consists of heat kernels, while the faster term con- 
sists of the hyperbolic part. Moreover they gave a more precise description of the behavior of the diffusive 
part, which is decomposed into four blocks decaying with different rates, and the conservative one. 
By using these refined estimates we were able to determine the asymptotic behavior of smooth solutions. 

The article is organized as follows. Next section deals with the modeling background relative to system 
(1). In the second section we recall some basic results about dissipative hyperbolic systems satisfying the 
Shizuta-Kawashima condition. In the subsequent section we show the global existence of small solutions 
by means of energy estimates. Finally, the last section is devoted to the study of the decay properties of 
small and smooth solutions to the quasilinear system (1). 

2 Modeling Background 

The formation of new blood vessels, called vasculogenesis, is a process lead by a chemotactic phenomenon, 
i.e. cells direct their movement according to certain chemicals in their environment. As a matter of fact, 
recent works [4, 18] have confirmed that endothelial cells in the process of vascular network formation ex- 
change signals by the release and absorption of Vascular Endothelial Growth Factor (VEGF-A). This growth 
factor can bind to specific receptors on the cell surface and induce chemotactic motion along its concen- 
tration gradient [13]. This communication by chemical signals determines how cells arrange and organize 
themselves. 

As shown in [29, 30], chemotaxis is decisive in many biological processes. For example, the formation of 
cells aggregations (amoebae, bacteria, etc.) occurs during the response of the populations to the change 
of the chemical concentrations in the environment. Moreover, also in multicellular organisms, chemotaxis 
of cells populations plays a crucial role throughout the life cycle: during embryonic development it is im- 
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portant in organizing cell positioning, for example during gastrulation [12] and patterning of the nervous 
system [31]; in the adult life, it directs immune cells migration to sites of inflammation [36], fibroblasts into 
wounded regions and during cancer growth it allows tumor cells to invade the surrounding environment 
[6] or stimulating new blood vessel growth [26] . 

This biological phenomenon can be described at different scales. For example, by considering the popula- 
tion density as a whole, it is possible to obtain macroscopic models of partial differential equations. One 
of the most celebrated model of this class is the parabolic one proposed by Patlak in 1953 [32] and subse- 
quently by Keller and Segel in 1970 [24]. 

However, the approach of PKS model is not always sufficiently precise to describe the biological phenom- 
ena [14]. As a matter of fact, diffusion leads to fast dissipation or explosive behaviors and prevents us to 
observe intermediate organized structures, like aggregations. This approach describe processes on a long 
time scale, while on a short time range one gets a suited description from models with finite characteristic 
speed. 

Kinetic transport equations describe quite well the movement of single organism. For example the "run and 
tumble" (the movement along straight lines, the sudden stop and the change of direction) can be modeled 
by a stochastic process called velocity-jump process [20, 35] . 

At an intermediate scale, the process can be described by means of hyperbolic differential equations. 
This class of models can be derived as a fluid limit of transport equations, but with a different scaling, 
namely the hydrodynamic scaling t — » et, x — » ex [5]. 

Hyperbolic models can also be obtained by phenomenological derivations, as done by Gamba et al. [15, 33] 
to describe the vasculogenesis process. In [15] the authors proposed a model including chemotaxis as a 
fundamental mechanisms for cell-to-cell communication in order to find key parameters in the complexity 
of the formation of vascular network. This biological process proceeds along three main stages: migration 
and early network formation, network remodeling and differentiation in tubular structures. The model 
proposed in [15] focused on the first stage of the process. 

The experimental results are encoded into a mathematical model starting from the assumptions that the 
cell population can be described by a continuous distribution of density p and velocity u, moreover it is 
also assumed the presence of a concentration (p of chemoattractant. The cell population in the early stages 
of its evolution can be modeled as a fluid of non-directly interacting particles and is accelerated by gradient 
of chemoattractant released by cells, which diffuse and degrade in finite time. 
From these assumptions follows the system 

d tP + V-(pu) = 0, 

< d t [pu) + V- {pu® u) + V_P(p) = -apu + ppV<p, (3) 
d t (p = D/S.(p + ap- b(p, 

where D, a, and b are, respectively, the diffusion coefficient, the rate of release and the inverse of the char- 
acteristic degradation time of chemoattractant, p measures the strength of cell response, and a the friction 
of the cells on the substrate. 

This system is derived in a classical way by Continuum Mechanics indeed the first equation describes mass 
conservation, the second one is a momentum balance with a chemotactic force and the last is a diffusion 
equation for the chemoattractant produced by endothelial cells and degrading in time; in particular, the 
convective term on the right hand side of the second equation allows describing cell migration. As indi- 
cated in [1] it is based on the assumptions that: (i) endothelial cells show persistence in their motion; (ii) 
endothelial cells communicate via the release and absorption of molecules of a soluble growth factor and 
this chemical factor can be reasonably identified with VEGF-A (Serini et al. [33]); (iii) the chemical factors, 
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released by cells, diffuse and degrade in time; (iv) endothelial cells neither duplicate nor die during the 
process; (v) cells are slowed down by friction due to the interaction with the fixed substratum; (vi) closely 
packed cells mechanically respond to avoid overcrowding. 

On the basis of experiments and theoretical insights, the authors in [33] showed that non-linear mechanics 
and chemotactic cellular dynamics fit into a model able to reproduce with great accuracy the formation of 
capillary networks in vitro. 

The model (3) is able to reproduce several experimentally observed facts, e.g. that the mean chord length is 
approximately independent on the initial cell density or that connected networks are formed only above a 
critical initial density as shown in [15, 33] . Moreover the authors provided a strong evidence that endothelial 
cells number and the range of activity of a chemoattractant factor regulate vascular network formation by 
flanking biological experiment, theoretical insights, and numerical simulations. 

3 System Properties 

Let us consider the Cauchy problem for the following hyperbolic-parabolic system 

d t p + V-v = 0, 

< d t v + V •(i><8>|) + VP(p) = -azJ + fipVcp, (4) 

d t cp - DAcp + ap - bcp, 
p(x,0) = poM, u{x,Q) - uq{x), (p{x,0) - (poix). (5) 

We made the assumption 

P'(p)>0 

to ensure the strictly hyperbolicity of system (4). 

Our aim is to prove that, under suitable assumptions, the Cauchy problem (4) -(5) admits a global smooth 
solution for small initial data. In particular, we are interested in solutions of the form (p, v, cp) = (p + p, v, cp + 
cp), where (p,0,cp) is a constant stationary solution to the problem with p > and cp = -|p, and (p, v,<p) is a 
perturbation. In this case, we can rewrite system (4) as follows: 

d t p + V • v = 0, 

* d t v + V-(v®-^Jj+\/P{p + p) = -av + p{p + p)\/(l), (6) 



d t <p = DA(f> + ap - bcp. 



Now, we show some properties of the hyperbolic-parabolic system (6). To this end, we rewrite it in the 
following compact form: 



d t U+ I d x ,MU+ U) = g(U+ U) + h(U+ U,V<p), 
7=1 



dt<p - DAcp + ap - bcp, 



where U = (p, v), U = (p, 0), fj(U+ W = \vj,^,...,^ + P{p + pl,...,^\, g{U+ U) = {0,-av) and h{U+ 



U,Vcp) = {0,p{p + p)Vcp). 
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3. 1 Strictly Entropy Dissipative Condition 

Let us consider the hyperbolic part of (6), that is 

d t p + V-v = 0, 

(7) 

d t v + \/-(v®-^Jj+VP{p + p) = -av. 

First of all, we want to prove that system (7) is endowed with an entropy function, that is a convex real 
function £ such that there exist related entropy-fluxes qj satisfying the following condition 

[fjYVS = q'j, 

where fj are the fluxes of the system (7), for j = \,...,n. 

To ensure the existence of entropy- fluxes qj, it is sufficient to prove that 

(fjYS" is symmetric, 

or equivalently that 8'' ' fj is symmetric, for each j = 1, n (see [3]). 

Once we have proved the existence of an entropy function for system (7) , the following additional equation 
for the entropy evolution can be written: 

d t {S{U) - VS{U) -U)+it d *j WW " V8{U)fj{U)) 

7=1 

= (V* (17] - VS{U)) ■ {g{U) - g(L0), 

where U is an equilibrium state for the system (7) (i.e. g(U) — 0). From this equation, we deduce that the 
integral of & (IT) - V<? (L/) • U is decreasing in time, if the term on the right-hand side is negative. 
Denoted by y the set of equilibrium states of the system (7), this property is encoded in the following defi- 
nition. 

Definition 3.1. An entropy & for the system (7) is dissipative at U, with U e y, if it satisfies the inequality 

{VS{U) - VS 0)) ■ (g((7) - g(0}) < 0, 

for any U in a neighborhood ofU. 

Clearly, an entropy function is dissipative, if it is dissipative at U, for each U £ y. 
Next, if 8 is a strictly convex function, we can introduce the entropy variable 

W:= V8{U), 

and the functions 

8*(W) := W-<$>{W)-g{<$>{W)), 
qj(W):= W-fjmW))-qjmW)), 

where <P = (Oi,0 2 ) := (V<£) _1 . Let us set now A = [S*)"(W), A, = f'A<t>[W))A and G{W) = g(<P(W)), then 

•* J 

we rewrite the system (7) in the entropy variable as 

n 

AodtW+Y, Ajd X] W=G{W). (8) 
7=1 
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Let us observe that the matrix Aq is symmetric positive definite and Aj is symmetric, for each j = 1, n. 
Now, we take an open subset of R" +1 and define 

y:={t/e^ : g({/) = 0}, 

r := V<?(y) = (»e Vif(t) : G(W) = 0} . 

Let us observe that the dissipative condition ensures the existence of a real positive matrix B = B(W,W) 
such that, for every W in a suitable neighborhood of W, 

G 2 [W) = -B{W,W){W 2 ~W 2 ). 

We refer to [17] for more details. 

Definition 3.2. The system (7) is strictly entropy dissipative, if there exists a real matrix B{W, W) e U n x U n , 
positive definite, such that 

G 2 (W) = -B(W,W)(W 2 -W 2 ), 
for every WEVSify) andW '= {Wl,W 2 )eT. 

In our case, we can consider for system (7) the canonical entropy function 



t 
P + 'P 

then 



1 V 2 _ fP + P P(T) 

S{p + p,v) = - - + {p + p)\ —y-dT, 

2 p + p Jo T 



S p {p + p,v) 



1 V 2 P(p+p) fP + P P{T) 



-p) r p+p pw, 

— + I 5-«T, 

P JO T z 



2{p + p) z p + p 

S v {p+p, v) = — V —. 

P + P 

It is easy to prove that system (7), endowed with the entropy <8, satisfies the strictly entropy dissipative con- 
dition. Moreover, according to [17], the definition of dissipative entropy is invariant for affine perturbation, 
so if we consider the function 

£(U) = g{U+U)-<g(U)-V£(U)-U, (9) 

we have that <§{U) is a quadratic function and still an entropy for (7). This system, endowed with S{U), 
satisfies the strictly entropy dissipative condition as well. 
Let us finally observe that, in our case, 

W := V<f ([/) = V@{JJ+U) - V<g(U), 

and 

y = {[/e^ : U={p,Q)}, 



r={WeV<#(<20 : W 2 = 0}. 
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3.2 The Shizuta-Kawashima Condition 

This section is devoted to prove that system (8) satisfies the Shizuta-Kawashima condition [34] that is: 

n 

(SK) every eigenvector of £ f'(.U)^j is not in the null space of g'{U), for every ifeR" - {0}. 
J=i ; 

Let us define A(U+U) = £ f'-{U+ U)tj, with £ e U n - {0}, so 
J'=i 1 

tn \ 




A(E7) = 



p'Go)^- o ... o 

\, P'(p)$„ o ... o / 

Now, we suppose that X is in the null space of g'(U), that is 

g'(U)X = <^ Xj = for j=2,...,ra+l. 

So if X is in the null space of g'(LT), then X = (Xi,0, . . . ,0) f . Since in this case we have 

XX = A(D)X <=^> XXi = 0, 

X cannot be an eigenvector of A(U). But this is absurd because X is an eigenvector of A{U), therefore X is 
not in the null space of g'{U). This proves that system (8) satisfies the Shizuta-Kawashima condition. 



4 The Global Existence of Smooth Solution 

In this section, by means of the entropy method, we aim to prove the global existence of smooth solution to 
the complete hyperbolic-parabolic system 



d t p + V-v = 0, 

d t v + V '•[■^p + PiP + P)) = -av + p{p + ~p)V(p, 
dt(p- D\(p + ap- b(p. 



(10) 



Let us recall that p, cp : R" x R + -^R + , u : R" x R+ — R", v.- [p + p)u, and P'{p + p) > 0. Moreover, U = (jo, 0, (p) 
is a constant stationary solution to the problem, with <p = f jo. 

As we have shown in previous sections, the hyperbolic part (7) of system (10) is endowed with the dissipative 
entropy (9) and it satisfies the strictly entropy dissipative condition. So, considering the entropy variable 
W = VS{U) and, setting <5(W) = (VJrW) and 



A {W) = mW))', 

aw) = gmw), 



Aj(W) = fjmW))A , 
H(W,V0) = h(®(W),V(f>), 



we write our system as: 



A d t W+ I Ajd Xi W=G(W) + H(W,V(b), 



d t (p = DAcp + aOi ( W) - b(p. 



(11) 
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Let us notice that the existence of a local solution to system (11) is ensured by classical argument. Indeed, 
(11) is a symmetric hyperbolic-parabolic system, therefore we know that, if initial data Wq, (po are in H s (U n ), 
with s > [nl2\ + 1, then there exists a local in time solution {W,cp) e C([0, T),H S (U")) x (C([0, T),H s (M n )) n 
L 2 {[0,T),H s+1 {R n ))) forsystem (11) (Theorem 2.9, [22]). 

Then, thanks to the continuation principle (see [27]), in order to prove that [W,<p) is global in time, it is 
sufficient to show uniformly (in time) estimates of the local solution. 
We state now our result: 

Theorem 4.1. Fix s > [nl2] + 1 =: sq. We consider the Cauchy problem associated to system (11), with small 
initial data Wq andcpo inH s (U n ). If\\Wo\\ H s, \\<po\\ H s and~p are sufficiently small, then there exists a unique 
solution {W,cp) of system (11), such that 

WeC{[0,oo),H s m n )), c/)eC{[0,oo),H s m' 1 ))nL 2 ([0,(x l ),H s+1 m n )), 

and for each t>0 



where C = C(\\W \\ H s,\<i> \ HS ,p). 

In order to prove this theorem, we firstly show the validity of some energy estimates for functions W and <p. 

Remark 4.2. This result holds for perturbation of small constant ( non null) states. Moreover in the isothermal 
case, i.e. P[p) = p, it holds also for perturbation of zero state in the one dimensional case. It can be proved 
using the result of [28] that ensures the existence of a dissipative entropy for 2x2 hyperbolic system, under 
suitable assumptions [11]. 

Remark 4.3. Concerning the constants, they all have been denoted by the letter c. Thus, c may stand for num- 
bers that are different from line to line of the text. Only when we intend to explicitly indicate the dependence 
ofc on some parameters, or to avoid confusions, we have used some other notations for the constants. 

4. 1 Energy Estimates for (p 

We consider the local (in time) solution {W,(p) e C([0, t),H s {K n )) x (C([0, t),H s m n )) n L 2 ([0, f) , H s+1 (R"))) , 
where te [0, T). 

First of all, we consider the parabolic equation 



Applying the spatial derivative of order a, with a = («i, a n ) and < |a| < s, and multiplying by d"(p, we 
get 





d t <p — DAcp + ap — b(p. 



(12) 



a f (i«3^) 2 ) + Df(a: 




a,+l 




.d^d^ + ad^pd^-Hd^f. 
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Then, integrating with respect to x and t , we have 

\^d a x <t» 2 dx +D |£ (d%d%Xj +1 ...dtyfyTdx 

~\f [d " ( l >o)2dx+ Ye\ Q f ( d "P^ dxdT 

+ lY~ b )f J ( d xW 2dxdT - 

Now, let us introduce the generic functional 

Nf(t):= sup ||W(t)|| 2 ,+ ( t \\W 2 {T)\\ 2 .dT+ [ WVWMWlr^dT, 
0<T<t "Jo "Jo " 



for / = 1.....S, and 



N${t):= sup ||W(t)||* 2 + f f ||W 2 (T)||2 2 dT. 

0<T<i JO 



Therefore, summing up the estimate of d"<p for a such that \a\ E [1,5] and assuming e sufficiently small, we 
obtain the s-order estimate for function cp: 

Mh' + cJ P^>fw dT + c j WVcpfjjs-rdTZ \\<f>of HS + cN*tf). (13) 

Remark 4.4. Let us notice that the previous inequality allows us to estimate the integral in time of ||V0|| RS 
simply using the functional N s {t), which involves the derivatives up to the order s ofW. 

Now, the function <p is a solution of the parabolic equation (12), therefore, using the Duhamel's formula, it 
could be written as 

(p(x, t) = e~ bt T p (x, t) * (p {x) + f e~ Ht ~ T) T p {x,t-T) * p(x, t)dr, 

Jo 

where T p is the heat kernel. Consequently, we obtain 

|V0(jc,f)|| ioo < £ [«~*'||3x,0o||. oo + sup ||p(T)|| LCO ( e- b[t - r \t-T)- llz dT 
i-\ L 11 " L o<T<t Jo 



<ce- bt {\\<t> \\ HS+1 + N s {t)). (14) 
4.2 The zero-order Energy Estimate 

Now, we want to estimate the L 2 -norm of the function W. To this end, let us rewrite the first two equations 
of system (10) in the form: 

n 

dtU+Y. dxjfjiU+ U) = g{U+ U) + h{U+ U,Vcp). (15) 

■/=! 

Multiplying the previous system by V<f (IT) = VS{U+U) - VS(U), we have 

n 

d t S{U) + £ d X] cjjlU) = V<# (U) ■ g{U+ U) + V<f ([/) • h{U+ U, V(/>), 
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where cjj {j = I,..., n) are the entropy-fluxes associated to the function £ . 

Let us observe that, thanks to definitions of the entropy S and the variable W = VS{U), there exist two 
constants So and c such that 

-\W\ 2 <<#([/) <c\W\ 2 , 
c 

for \W\ < 8q. Moreover, as proved in Section 3.1, the system (15) satisfies the strictly dissipative condition, 
therefore there exists a constant c such that 

-(W-G(W)) > c\W 2 \ 2 . 

Let us integrate the previous system, with respect to space variable x, so we get: 

Y t j £(U)dx = j Vi(U)-g{U+U)dx + j VS(U)-h(U+U,V(p)dx, 

which yields 

\\W(t)f L 2 + cj \\W 2 (T)\\ 2 L2 dT<\\Wo\\ 2 L2 + cJ o J Vg{U+U)-h{U+U,V(p)dxdT, 
for all \W{x, t)\<8 where (x, t) e U" x (0, T). 

Now, let us observe that, thanks to the definition of function h{U +U, Vcp), the last integral can be estimate 
as follows 



J J WS{U)-HU+U,V(l))dxdr=J J W 2 - y.(p + ~p)V(pdxdT 

< [ n\\W 2 {T)\\ L 2\\ P {T)\\ L ~\\V(P{T)\\ L zdT 

Jo 

+ [ A*P II W2 (T) || x.2 || V0(T) || i2 rfT 

Jo 

<c sup ||p(t)||l« [ {\\W 2 (T)\\ 2 L 2 + \W(l>(T)\\ 2 L 2)dT 

T£(0,f) JO 

+ Cp J ' [ || W 2 (T) || 2 L2 + || Vcp (T) || 1 2 ) dT 

< cTVi (?) [ 1 0o || % + Nf ( t) ] + cp [ || 00 1| 2 L 2 + N 2 ( t) ] , 

where, in the last inequality, we used the energy estimate of the function (p. 
In conclusion, the zero order estimate of function W is given by 

II W(0 II 2 L z + f II W 2 (t) || 2 L2 dr < N 2 (0) + CC|| (po || ( 



+ i 



C(p)N 2 {t) + CJVf (f) + C(||0o|| £2 -P)- 



4.3 The 5-order Energy Estimate for function W 

To prove energy estimates of the function d s x W, it is necessary to use some inequalities based on the Sobolev 
embedding theorem. Here we just state the following lemma, whose proof can be found in [37] . 

Lemma 4.5. We take s, s\ ands 2 three non-negative integers and sq := [nl2] + \. Then 

(i) if S3 = min{5i,s 2 ,si + s 2 - s } > 0, then H Sl H S2 c H S3 

(the inclusion symbol c denotes the continuous embedding); 
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(ii) if s > sq and A', U e i/ 5-1 , then for all multi-indices a with < \a\ < s, the commutator [d®,A]U := 

d«{AU) - Ad a x U e L 2 and 

|| [d?,^][/|| L2 <c s ||^V-i ll^llfli«rt-i; 

(iii) ifs > sq, V e H s with values in Q, and A e C S (Q) with A{0) = 0, then A{V{-)) £ H s and 

\\A(V(-))\\ H s < C s \A\ s \\V\\ H s a + \W\\^}). 

Here, C s is a constant depending only on s and n, and 

\A\ S := sup \dfjA{U)\. 

LTen,l<|a:|<5 

Nowwe estimate the L 2 -norm of the s-order derivative of the local function W. To this end, we consider the 
system 

n 

d t W+ Y,Ajd Xj W=A - 1 G{W) + A - 1 mW,V(l>), (16) 

where Aj :- A^Aj. Applying the derivative d x , where 1 < \ a\ < s, we get 

d a x d t W+£ i A J d%d Xj W= A-'d a x G+\d a x ,A-']G+A-'d a x H 

Hd a x ,A- l ]H + Y. n j=l [Aj,d«]d Xj W, 

where [a, b]c := a{bc) - b{ac). 

If we multiply this equation by {d" W) t Ao, we have 

n 

{dlwYAod^dtW+Y, id" W) 1 Ajd a x d Xj W = {d"W) t d x G+ (d a x W) * Ao [d x ,A Q ~ 1 ]G 
/=i 

n 

+ (d a x W) 1 d a x H + (d a x W) * Ao [d a x , Aq 1 ] H + £ {d%* W) ' A [A ) , d a x ] d Xj W. (17) 

i=i 

Thanks to the symmetry of Ao and Aj , we deduce the following equalities 

{d a x WY A Q d a x d t W ^d t [{d a x WY A Q d a x W)- ^{d a x wY {d t A Q d a x W) , 

£ [d a x wY Ajd a x d Xj w =\ £ d Xj {{_d a x W) t A J d a x w) -If; [d a x wYd Xj Ajd a x w. 

Let us observe that, thanks to the strictly dissipative condition, there exists a positive definite matrix B such 
that 

[d a x WYd a x G= -{d a x W z Yd a x {BW z ) = -{d a x W 2 Y B{d a x W 2 ) + [d a x W 2 Y[B,d a x ] W 2 . 
Substituting these equalities in (17) and integrating with respect to the space variable, we obtain 

\~ahf ( d " w ^ A o d " wdx + J id a x W 2 YBd a x W 2 dx^ j {d a x W 2 Y[B,d^]W 2 dx 

+ [ {d a x WY[Ao[d a x ,A- Q l ]G+j^Ao[Aj,dl\d Xj W + Ao[d a x ,AQ l ]H\dx 
J \ j=i ) 
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- j (d a x W) T d t A +Y,dxjAj d a x Wdx + j {d a x WYd a x 
2 J \ j=i J J 



Hdx. 



(18) 



Let us analyze these integrals separately. Some of them can estimate by classical arguments, following the 
approach of [ 1 7, 37] . 

Above all Aq and B are positive definite, so there exist two constants such that 



{d%WY A Q d a x W > c\d"W\ 2 , 



(d% W 2 ) l Bd a x W 2 > c\d% W 2 \ 2 . 



(19) 



These inequalities allow us to estimate the left-hand side of (18). 

Now, we can estimate the time integral of the first term on the right-hand side of (18) using Lemma (4.5). 
Indeed, thanks to the condition (ii) and the regularity of functions, we have 



f*J \ { d« x W 2 YlB,d° x ]W 2 \ dxdr < cj \\b'\\ hs -i [\\d x W 2 \\ 2 L2 + llW^ll^ai.jj dr. 



(20) 



Then, we consider the second and the third integral on the right-hand side. We know that d" W e H s ~ w 
and Aq e H s , so 53 := min{s,5- |a|,2s- |a| - sq] is a positive constant. Therefore, using the condition (i) of 
Lemma (4.5), we obtain that A d x W e L 2 and ||A)d"W|| L2 < c||a£ W|| HsHa | IIAillfF- Then, using again con- 
dition (ii) of the same lemma and the regularity of functions, we deduce that [d x , Aq]G e L 2 and [d x ,Aj]d x , W e 
L 2 , for each j - 1, n. As a consequence of these observations, we can calculate 

J J ^"WY A Q [d a x ,A^ l ]G\dxdT 



< c 



(21) 



and 



17 



dxdr 



< c f \\A \\ H s £ \\a'j\\ (\\d a x W\\ 2 H ^ a] + ||VW|| 2 H|[IH1 ) dT. 



(22) 



Moreover, in the same way, we get 

|(a;w)%[a;,V]H| dx<c\\d^w\\ HS - M \\a \\ H s ||(V)'L-i u^iim-i-i 



/I 



Let us point out that classical arguments are not a sucessful strategy to estimate the r.h.s. of the previous 
inequality. Since the last term causes the failure of standard approaches, our aim is to show an effective 
technique to estimate it. 

So, we focus our attention on this term and we get the following inequality: 

jH{U+ U, V0) J i < /i [ II P II L « II V 4> II HM-i + II P II ffM-i II V( t> II l» + P II V< P II fl1«Hi ] • 
Now, substituting this estimate in the previous one and integrating with respect to t , we get 

J J \{_d%W) f AQ[d%,A^]H\dxdr 



Jr. f 
[\\ da x w fw-\<*\ + P<pf H \«\-i)dT 
v ' 

+ SUp ||p(T)|| H | ahl ( (||32W||^ J _| (t | + ||V^l|^ 0o )rfT 
0<T<f J0 v ; 
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+ p/ o f (lKw|| 2 HsH£[| + ||V0||^ a| _ 1 )] 
< c A [ N s ( t) ( AT S 2 ( t) + || c/) 1| z m + N 2 a ( f )) + N a - 1 (TV 2 ( r) + || O || ^ + N s ( «) 

+ p(wf(r)+|| ( /)o|| 2 H , tI i + ^(^)] 
<c^(7V s 3 U) + c(||0o|| HI )JV s (f) + c(p)N 2 (t) + c(||0o|| H i«i »P))» (23) 
where 71:= sup \\A \\ H s || (^" 1 )'|| ffs _i. 

0<T<t 

Next, the last integral of inequality (18) can be studied in the following way: 

<MO3jW 2 OLa(Oplli«-OV0||fl« + Opllfl«l|V^|| L .) 

+Afplia^w 2 || z2 ||V0iiH«, 

which, integrating with respect to the time variable, yields 

f f (d?W)^HdxdT<p SU P Upllico f f (||^W 2 || 2 2 + ||V0||^)dT 
JO J 0<r<t JO 

+ p SUp Upllfl* ('(||a?W 2 || 2 +||V(/)|| 2 e„)dT 

0<T<t JO 

< c7V s 3 (f) + c{\\(p \\ H s,p)N s {t) + c{p)N 2 a {t) + c{\(f) \ Ha ,p). (24) 

Remark 4.6. Let us point out that, in order to estimate the second integral of (24), it is not useful to consider 
sup || V(/>|| loo . Indeed, it is impossible to estimate this term by the functional N s , since, as deduced by (13), 

0<T<t 

the order of the functional should be increased up to s+\. While, as noticed in Remark (4.4), we can control 
the time integral of ||V0|| iOO by N s , without increasing the higher order derivative. 

Now, we examine the remaining term of (18). Using (16) and the definition of A;, we can write 
dtAo+f^d Xi A ] = -A^\f^A ] d Xj W + A' {A^G) + A' {A^H) + £ A'jd Xj W 



= I>o Ufa w) + A' {A^G) + A' [A^H) . 

7 = 1 



From this equality, recalling that G = (0, -BW 2 ) 1 , we deduce 



/ 



dtAo+Y.dxjAj 
;'=i 



d°tW 



I n 



dx<cA 



E PxjWj 2 + \\W 2 \\ L 2 
Vj=i" 



\ d " w \\h 



+ c^(||(p + p)V0|| L2 )||^W|| 2 2 
<cA( || W 2 \\ L 2 \\d a x W\\ 2 L2 + \\VW\\ L 2 \\d a x W\\ 2 L2 
+ ||(p + p)V0|| i2 |Kw|| 2 2 ), 



(25) 



where A:= sup < \\Ao\\ H * I U' + |K|| HS -i || A,, 1 ! „ s (1+ \\B\\ H >) 

0<T<t [ j = l 11 1 '' H 

Let us analyze these terms separately. First of all, we have 



rt r t 

\\W 2 \\ L 2\\d a x W\\ 2 L2 dT< sup \\d a x W{T)\\ L2 [\\d a x W\\ 2 L2 + \\W 2 f 2 )dT 

JO 0<T<f Jo v ; 



A HYPERBOLIC-PARABOLIC MODEL FOR VASCULOGENESIS 



15 



and 

[ l \\\/W\\ L 2\\d«W\\ 2 L2 dT < N a (t){Nf{t) + Nl^it)). 
Jo 

Now, we are interested in studying the last term of the inequality (25). Let us observe that 
f Up + pW<p\\ L 2\\d%W\\ L 2dT< sup \\p\\ L 2 [ [\\V(p\\*+\\d%W\\ 2 L2 )dT 

JO 0<r<t Jo v ' 

+ SUp ||p(T)|| LOO ['(WVcpllz + WdtWfJdT 
0<r<t Jo v ' 

+ P^(PH 2 L^Kw\\ 2 L 2)dT 
<c[\\(p \\ H s)N s {t) + c(p)N 2 {t) + N^t) + c{\\cp \\ L2 ,p). (26) 

Finally, integrating the equation (18) with respect to the time variable and substituting in it inequalities (19), 
(20), (21), (22), (23), (24) and (26), we deduce 

|| 5 a x W{ t) || \ x + f || d% W 2 (T) || \ 2 dT < C(|| W 1| H s , || 00 1| H s , p) + C(]| 00 1| H s ,p)M s (t)N s {t) 
JO 



where 



M s {t):= sup 

0<T<f 



+ c(p)M s (t)N^t) + cM s (f)AT s 3 (f), 



ii^oiih. + II (V)'llff. + ii^oiih* E \\ A 'j\\ ws , + ii^oiiff* IK^'llif-i 

y=l" 



(27) 



+ ||b'|| H5 -i + IKL-i IK^o^ll^ a+\\B\\ H °) 
+ (i+ iiwii hs -i + iiwn^) 2 wbw 2 ^ + 1 + ii w\\ H s-i + www 5 -^ 

+ (tW0)-Aj(W)\\ H s-i] 

Therefore, summing up for 1 < \a\ < 5, we deduce the following s-order estimate of function W 

II W{ t) || % s + [ || W 2 (t) || % s dr < C( || W II , || O || H s , p) + C( || O || w , p) M s ( t)N s {i) 
Jo 

+ C[p)M s [t)N 2 {t) + CM s {t)N^t). 
4.4 Proof of the Global Existence Theorem 

Now, we are finally able to prove Theorem (4.1), showing the existence of a global smooth solution for sys- 
tem (11). 

Proof. Let us recall the definition of the functionals 

N 2 At):= sup ||W(t)|| 2 „,+ f \\W 2 {T)\\ 2 ffl dT+ [ WVWMW 2 dr, for /= l,...,s, 
0<T<t "Jo " Jo " 

N^{t):= sup ||W(t)||2 + [' \\W 2 {T)\\ 2 2 dT, 
0<r<t Jo 



A HYPERBOLIC-PARABOLIC MODEL FOR VASCULOGENESIS 



16 



and the energy estimates, obtained in previous sections, 

IIW(fllli2+ f IIW2(T)||2 2 dT<A^(O) + C(||0o|| i2 )JVi(fl + C(p)iVf(f) 

JO 

+ CAf 1 3 (0 + C(||^ || L 2,p), (28) 

and, for each s > 1, 

II W( t) || % + f II W 2 (t) II %, dT < C ( || W 1| H s, || O || , p) + C( || 0o || ffs , p) M s ( f ) N s ( r ) 
Jo 

+ C(p)M s (t)N 2 it) + CM s (t)N*(t). (29) 
Therefore, to obtain an estimate of the functional N 2 {t), we have to study also the term 

J WVWWf^dT, for J = 1, 
To this end, we rewrite the first equation of system (1 1), in the following way 

n 

d t W + Y, A i W)dx] W = Aq 1 ( W) G{ W) + Aq 1 ( W) H(W, V0) + L{W, d x W) , 

where L := X (0) - Aj {W)) d x W. Applying the Fourier transform with respect to x, we obtain 

7=1 ' 

n 

d t W + i £ Zj Aj (0) W = A~ 1 G + A' 1 H + L, (30) 

Let us recall that, in Section (3.2), we proved that the first equation of system (1 1) without the term H{ W, V0) 
satisfies the condition (SK) . As shown by Shizuta and Kawashima [34] , this means that there exist a constant 
c> and a skew- symmetric real matrix K - K{^) e C°°(S" _1 ) satisfying K{-£) - -K{^) and 

i + [K({)A(0Y] + |f|dlag(0,/») > c\S\I n+ i, (31) 

for every £ e S" -1 , where S" -1 is the unit sphere in U" and 

A(0:=YAj(.0)Sj, {cR n \{0}. 

Now, if we multiply the system (30) by - iW l K, then we have 

n 

- i W' Kd t W + W f K Y Z j Aj (0) W = - i W'K(A~ 1 G + A~ 1 H + 1) . 

Substituting inequality (31) and 

2ImW t K(Aj r G+Aj r H + L) < c\{\\W\ 2 + C\^~\\ A^~G\ 2 + \A Q T H\ 2 + \L\ 2 ) 
in the previous system, we obtain 

-id t [W'KW] + c\{\\W\ 2 < 2|{ || W 2 \ 2 + Clff 1 (| ^Gl 2 + \Aq 1 H\ 2 + \L\ 2 ). 
Let us multiply this last inequality by |£| 2fc ~\ with k > 1, and integrate over U" x [0, f], so we obtain 

c£j \{\ 2k \W\ 2 d$dT<2f* J \t\ 2k \W 2 \ 2 didT 
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\tf k - l \W((,t)\ 2 df+cj Ifl^lWbl 2 ^ 



+ c 



Then, since 2|^| < l + k r | andfc>l, we deduce that 



2 ([A" 1 G| 2 + \A~ l H\ 2 + \'Lf)d(dT. 



f E \\d%WM\\ 2 L2 dT<c ft E ||<5?W 2 (T)|| 2 2 ]dT+||mf)ll 2 Ht + IIWoll 2 H , 

•'O \a\=k [ J \\a\=k ) 

+ f E (KcVa|& + ||a;cVfl)lM, 

J |a|=fc-l 



+ f E 11^ 

J |a|=fc-l 

which, summing over all a such that \a\ £ [1,5], yields 

| Q f II V WCt)!! 2 ^! dT <c f ^ || WzWH^ + || W{t)\\% + II W f H s 



+ j o (HVgII^-x + IIV^IIh-i+II^Ih-i)^ • 

Now, let us recall that £ = Z (-Ay (0) - Aj (W)) <3 X) W, so using condition (i) of Lemma (4.5), we get 



/'=i 



liP-i <cE II A/ (o) - Aj{W) || flS -i a Iy w 



then 



filing 



_idT < cM s {t)N z (t), 



where M s (f) is defined by (27). 

Using again Lemma (4.5), we deduce that 

|| Aq 1 G|| H _! < || 1 (0) G|| Hf _! + || [Aq 1 ( MO - 1 CO)] G|| h _i 

< c (1 + || Aq 1 (WO - Aq 1 (0) || || G|| H _i 

< c(l + ||W|| HS -i + IIWH^JIIBII^-i M^-i , 

which yields 

I J Al 1 G || jp_, < cM s ( t) N 2 S _ l it). 
Jo 

Proceeding in the same way, we get 

|| V zr || H _! £ cU + UWUh-i + IIWB^^iiHilfp-i , 

so, we have 

J || Aq 1 I \,. x dT<M s {t)J^ (|| p || LOO || V0 1| ffs _! + || p || HS .! || V0 1| ioo + p || V0|| HS -i f dT 
<M s {f)\ sup ||p||^oo f ||V(/j|| 2 idr 

L 0<T<t JO 
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|| V0 || 2 LOO dr + p 2 I || V0 1| HS _i 
o Jo 

+ 2 sup ||p|| i00 sup IpI^-i f (||V0||^ s -i + ||V0||^ oo )rfT 

0<T<f 0<T<t JO V ' 

+ 2p sup \\p\\ LOO I ||V0||^ s -idT 

0<r<f Jo 

+ 2p sup \\p\\ H ^ f (|| V0||^ +||V0||^-i)£/t] 
0<r<f Jo v ; J 

M s ( f ) [ N 2 ( t) ( || O || ^ . + JVf_ ! ( fl ) + N 2 ( fl ( || O 1| ^ + N 2 ( fl) 

+ P 2 ( || 0o || % + N 2 (t)) + 2N 2 (fl (fl O || 2 HS + N 2 ( t )) 

+ 27V S 2 ( t ) ( || 0o || %s + N 2 ( t )) + 2p7V s (fl (fl 0o || ^ s + N 2 ( t )) 

+ 2pN s ( f ) ( || 0o || ^ s + N 2 ( t )) + 2p AT S _! ( f ) ( 1 0o || 2 W + N 2 ( fl) ] 



<c(||0o|| ffs ,p)M s (f)Af s 2 (f) + M s (f)A^(D + c(||0o|| ffs ,p,M s (f)) 
+ c(||0o|| ffs ,p)M s {t)N s {t) + c(p)M s (t)N*[t). 



Consequently, as long as M s ( f ) < C, we obtain 

J" II V W || ^ dr < cN* (0) + c( || O || ffs , p) M s ( f ) N s ( fl + c( || 0o || Hs , p) AT, ( fl JVf ( fl 

+ c{-p)M s (t)N]{t) + cM s {t)N*{t) + c(||0o|U ,P>M s (t)). 

Combining the previous inequality with (28), (29), we get the estimate 

N 2 (t) < CN 2 {0) + C(||0o || H s ,p,M s (fl) + C(||0o \\ H s ,p,M s (fl)AT s (fl 
+ C(p,M s (t))N 2 (t) + C(p,M s (fl)N s 3 (fl + C(M S ( fl) AT s 4 (f). 

In conclusion, choosing small initial data and small constant state, from the previous inequality we deduce 
the theorem, by classical arguments. □ 



5 Asymptotic Behavior 

In this section we study the time decay properties of the global smooth solution to system (10), proceeding 
along the lines of [2] . Thanks to the decomposition of the Green function of the linearized problem, we aim 
to obtain the H s and L°° decay estimates of the solution for the considered model. 
To this end we rewrite system (10) in the Conservative-Dissipative form as 



d t {U+U)+ I d x .fj(U+U) = g(U) + MU+U,V<p), 
7=1 



df = DA0 + ap- b(p, 



where 





t p 1 








u= 


V 


, u= 


{ j 


, fj(U+U) = 




I Vp' (pi J 







VP^Vj 



'iv + P(p+p) c 



p+p 



(32) 
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-av 



h[U+ U,V(f>) 





P+P 
v/f"(p) 



V</> 



Defined fj(U) = fjiU+U)- fj (LO and £ = 



AP'(p) 



, the system can be rewritten in the following way 



(33) 



(34) 



dtU+Y, d Xj [f'j(U)u) = g(U) + £ 5 X . (/^I7)tf-/j(t/)) + h(U+ U,V<f>), 

7 = 1 7=1 

and its solution is given by 

U{t) =T h {t) *U +f^f o d Xj r h {t-T) * [7^(I7)t/(T)-7 J -(t/(T))] dT 

+ [ T h {t-T)*h{U+U,V(p)dT, 
Jo 

where r'' denotes the Green function of the linearized system 

d t U+Y J j' ] {jj)d Xj U = g{U). 

7 = 1 

Let us briefly recall the results on the Green Kernel of multidimensional dissipative hyperbolic systems 
obtained by Bianchini et al. in [2]. In their work the authors analyzed the behavior of the function T (x, t) 
for linearized problems. It has been decomposed into two main terms: the diffusive one consisting of heat 
kernel and a faster term consisting of the hyperbolic part. In general, the form of the Green function is not 
explicit, but it is possible to deal with its Fourier transform. The separation of the Green Kernel into various 
parts is done at the level of a solution operator T (t) acting on L 1 (IR") n L 2 {U n ). 

They deeply described the behavior of the diffusive part, which is decomposed in four blocks, decaying with 
different rates. They showed that solutions have canonical projections on two different components: the 
conservative part and the dissipative part. The first one, which formally corresponds to the conservative 
part of equations, decays in time like the heat kernel, since it corresponds to the diffusive part of the Green 
function. On the other side, the dissipative part is strongly influenced by the dissipation and decays at a 
rate t~z faster than the conservative one. 

They considered the Cauchy problem for the linear system in the conservative-dissipative form 

n 

dtW+Y^ Ajd Xj w-Bw, 

7 = 1 

and they showed that it is possible to decompose the solution as 

w(t) = Y h {t) * w - K{t)wo + JC{t)wa, 

for any function wq e L 1 (IR") n L 2 (U n ), where K{t] is the diffusive part and J£{i) is the trasport dissipative 
one. 

Moreover for any multi index /3 and for every p e [2, +oo] the following estimates hold: 



IID^Wu/olliE <Ce- cf ||D^wo|| i2 , 
WLqDPkWwoWlp <C(|)8Dmin{l,r f ( 1 "p)"^}||LoM'ollLi 
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+ C(|/J|)min{l,f 2 1 pI 2 2 H|L_u^|| L i, 



HL-D^^COm/oIIip <C(|/3|)min{l,r "i 1 ') * 2 }\\L w \\ L i 

+ C(|/3|)min{l,f 2 1 1 pi 1 2 }||L_w || L i, 
where Lo = Ui , 0] and L_ = [0, .fe] are the projectors on the null space and on the negative definite part of B. 

5.1 H s Estimates of the Solution 

This section is devoted to study the decay rates of solution to the system (32) in the H s -norm. 
We define 

E s :=max.{\\Uo\\ L i,\\Uo\\ H s}> D J :=max{||<j&o|| L ii WIh*}' 
and the general functional 

S%:= sup {max{l,T a }||w(T)|| H 4. 

0<T<t 

Then, we shall prove the following theorem 

Theorem 5.1. Let {U, <p) be a global solution to problem (32), with initial conditions 

U(x,0) = Uq(x), (p(x,0) = (j) {x), 

with 

UnE H s+1 m")nl}m"). (hne H s+1 (R n )nl}(R n ). fors>\ 



Lf eH s+1 (lR")nL 1 (lR"), (f>o e H s+1 [R n ) n L 1 (M n ), /ors>[-^]+l. 
Then the following decay estimates hold: 



|| U ( t) || h. < minU , r t } C(E S+1 , D s+1 , p) , 
110(f) || H s + i < min{l, ri}C(E t+1 + D,+i,p). 

Proof. First we consider the parabolic equation 

d t <p = DA<p + au - b(p, 

and, using the Duhamel's formula, we can write the solution as 

(p(x, t) = [e~ bt Y p {t) * O ) (x) + f e~ b{t ~ r) Y p {t-T) * ap(T)di, 

Jo 

where 

e iDf 

T p (x, t) : 



{AnDt) nl2 ' 
Let us start with the H s+1 estimate: 



\\(pm H ^ <ce- bt \\<p Q \\ H s+i +c [ e- b[t - T) \\ap{T)\\ L 2dT 

Jo 

+ c[ e~ blt ~ T) {t-Tr^\\ap{T)\\ H sdj 
Jo 
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<ce- bt \\(p \\ H s + i +cs|cfl f e- &(t - T) (f-T)-2min{l,T-7}dT 

JO 



apt l 
+ cS^(t) J e" 6(f " T) min{l,T"4}dT. 



So we obtain the following H estimate for <p 



\\c/){t)\\ H s +1 <c(e- &t ||0 o ||^+i+min{l,|r-ir7}S*(r) + min{l,r7}S^(f)), 



(t)<C(e~* t max{l,^}||^ || fl w 1 +S^(t)). (35) 



which yields 

Let us notice that from the previous inequality the decay rate of the function (p in H s+1 is the same rate of 
the function U in H s . 

Proceeding in a similar way, it is possible to get the following estimate for the function (p in the space L 1 : 

||0(r)|| L i<e- fc '||0ollLi + c sup llpWIlii, (36) 

TE(0,f) 

where, thanks to the mass conservation, sup ||p(T)|| L i = ||poll{i. 

TE(0,f) 

Now we focus on the estimate of function U. Let us observe that fj{U) - f AU)U - U 2 rj{U) (where the 
product should be intended as the tensor product), therefore, using (34) and the definition of E s , we obtain 

II U{ t) || H s < c min{ 1 , f~ * } II Uo II L x + ce~ c f \\ U 1| H s 



nt j n 

+ c i min{l,(f-T)"4-2}^ \\U 2 (T)rj{U(j))\\ Ll dT 
Jo 7= i 

[' e -c(t-i) y ||a x .(Lf 2 (T)r,-(Lf)(T))|] rfr 



+ c 

+ [ \\T h (t-T)*MU+U,V(l))(.T)\\ H sdT. (37) 
Jo 



At this stage we want to estimate the right hand side of this inequality. 
Let us start studying the first integral in (37), as follows 



( min{l,(f-T)-f-2}£ \\U 2 (.T)rj(mT))\\ Ll dT 
Jo ]=1 

< ( f min{l,(f-T)-7-2}||[/(T)|| 2 2 £||ry(t/(T))0 B » flD1 <So) dT 

JO j = l 

< c(S^(f)) 2 J min{l,(f-T)"^"2}min{l,T"t}d T . 



Then from Lemma 5.2 of [2], we deduce 



c f min{l,(f-Tr?~^}||f/ 2 (T)r/([/)LidT 
Jo 

< cj min{l,(f-T)"7-i}min{l,T"?}(S^(f)) 2 

<cmin{l, f" v }(Sj(f)) 2 , (38) 
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where v = min{^ + |, §, §«- \}. 

In order to estimate the next term in (37), we use Lemma 5.3 of [2] which yields 

£ \\d Xj (U 2 rj(U))\\ < sup c(6oA\u\\ H sA\rj\\ C s m[M<Sa) )\\U\\ L ~£. l a *j u l„. 
;=l" j=l,...,n y=l" 

<c||I^||fl.||£?|| flW i. (39) 

Then we have 

f'e-cw £ ||a x/ ([/ 2 r,(C7))(T)|| <c f e"*^ ||U(t)|Ih* II 17(t) || h «.i dt 

JO j =1 » »H" Jo 

<cS*(t)E s+ i j e" c(f " T) min{l,T"4}d T 

J 

<cmin{l,f-7}sj(f)£ s+ i. (40) 

In the last inequalities, we have used Lemma 5.2 of [2] and the estimate of Theorem (4.1) to controll the 
norm of the function U in H s . 

Finally we estimate the last integral of (37) in the following way 

[ \\Y h {t-T)*h{U+U,V(p){T)\\ H sdT<l \\JC(t-T)h(U+U,V<p){T)\\ H sdT 
Jo Jo 

+ f \\K(t-r)HU+U,V<l>KT)\\H'dr. 
Jo 

For the first term, we have: 

f \\jrtt-T)HU+U,V<l>)(T)\\HsdT< [' Ce- Clt - T) \\V(!>{T)\\ H s(p+\\p(T)\\ H s)dT 

Jo Jo 

<7iS^ (r) J ce~ c[t ~ T) min{l,T~7}dT 

;cmin{l,r5}ps| x (0 



+ S^(f)S r 4 r (r) / ce" clt_Tj min{l,T"2}dT 



+ cmin{l,r2}S^(r)S^(r). (41) 

In order to complete our estimate, we need to study the contribution of the diffusive part of the hyperbolic 
Green function. Since we are interested in the slowest decay estimate of the solution U, we focus on the first 
component: 

f \\K(t-T)h{U+U,V(p){T)\\ H sdT 
Jo 

>t n 



C ii 

I L \Ki j+ i{t-T)d Xj (p{p + ~p){T) \\ dT 
Jo 

cf min{l,(f-T)"'J- 1 }7J||0(T)|| L idT 
Jo 



Thanks to (36), we deduce that 



cf min{l,(f-T) " l }p\\(p(T)\\ L idT<cmm{l,t " Hpll^ollz 1 + c\\p \\^pt (43) 
Jo 
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In conclusion, substituting (38), (40), (41), (43), in (37), we have 

\\U{t)\\ H > <c[mm{lJ~hE s + mm{l,t~hslE s+1 +mm{l,t~ v }(sfj(t)f 
+ min{l, r i }pS ^ ( t) + min{ 1, r f } S £ ( S * ( f ) 
+ min{l, r2}sj a (flsj(t) + /Zmin{l, r? _1 }p||^o|| £ i + llpoOzipf 

So we obtain 

S I ( t ) < c (j?, + S I ( r)£ s+ 1 + (s| ( fl ) 2 + ps| x ( + pD s + s| v (?) S lit)] . 
Now, we substitute inequality (35) in the previous one, obtaining, for t > 8 > 0, 

Slit) <C{\ + Sl{t) + [Sl{t)) 2 ), 

where C = C(f! s ,D s+ i,p). 

From this inequality we deduce that, if the initial data and the perturbation p are sufficiently small, then we 
have 

\\Uit)\\Hs<Cmin{l,rh, 
||0(t)|| HS+1 <Cmin{l,r?}. 

5.2 L°° Estimates of the Solution 

We now estimate the L°°-norm of solutions to the system (32). As done before, we define the functional 

Rlit):= sup {max{l,T a }||w(T)|| i( »}, 

0<T<f 

and 

E s := max{ || U \\ L i , \\ U \\ H s} , D s := max { || O || £ i - || </>o || H , } ■ 
We want to prove the following theorem 

Theorem 5.2. Let iU, 0) be a global solution to system (32), with initial conditions 

Uix.O) = U ix), 0(x,O) = (poix), 

with 

U eH s+1 [« n )nL 1 (.« n ), O G-ff s+1 (K")nL 1 (R"), fors= [-] +2. 



□ 



Then the following decay estimates hold: 

\\Uit)\\ L oo<mm{l,rhciE s ,D s+l ,p), ||0(f)|| L °° <min{l, r? }CiE s ,D s+1 ,p). 
Proof. Proceeding as done before, we obtain L°° estimates for and V0. First of all we show that 

|| 0( t) || L oo < ce~ bt || 0o II l« + c f ' e- fc(f - T) || ap (t) || l °° 

Jo 
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which yields 



<ce- ht \\((>oh°° + cR^tt) f e- &(t - T) min{l,T-f}dT, 

JO 

\\(p[t)\\L<*> < c(e" &t ||(^o || L ~ + min{l, H }flj(t)) . 



In a similar way, we get 

|| Vcp ( f ) || LOO < c [e~ bt || V0 O II + min{ 1 , 1 1 - 1 1" * }R% it )) . 

This means that 

Rl<C{D s + Rl{t)), (44) 
fl| a <C(D,+i + flJ(t)). (45) 
Now let us consider the solution of our system written in the form (33), i.e. 

U(t) =r'\t) *u +jrJ o d Xj r h {t-T) * \ f](U)U{T)-fj{U{T))\ dr 

+ [ T h (t-r)*HU+U,V(p)(r)dr. 
Jo 

Thanks to the decomposition of the Green function, we estimate the L°°-norm of U in the following way 
|| U( t) || ico < cminfl, t~ % } II U || £ i + ce~ ct || U 1| H * 



J-t j n 

min{l,(f-T)"f-2} £ \\U 2 rj{U){T)\\ ! dj 

+ [ e - c{t - r) j^\d x AU 2 rj{U){T))\ dT 



10 

The third term in the r.h.s. of (46) is estimated as: 



+ f \\T h (t-T)MU+U,V(!)){T)\\ L °odT. (46) 
Jo 



I min{l,(f-T)"f"2} £ ||[/ 2 rj([/)(T)|| Ll dT< cmin{l,r^}£f 
Jo ;=1 



While the next term in (46) can be estimated as 



£ f e -cf-T) II a [C/ 2 cc ^ Ct)) II dT<cR^t)E s+1 f e- ctf - T) min{l,T-f}dT 
j = ]_J0 " Jo 



<cmin{l,r4}i?*( ? )£ s+1 , 



where we have used the inequality (39) . 

Proceeding in a similar way, we control the last term in (46) as follows 



[ \\T h {t-T)*h{U+U,V<p){T)\\ L c°dT<{ \\K{t-T)h{U+U,V<p){j)\\ L <x>dT 
Jo Jo 

+ f \\Jr(t-T)h(U+U,V(p)M\\ L oodT. 
Jo 



A HYPERBOLIC-PARABOLIC MODEL FOR VASCULOGENESIS 



25 



Let us start from the second integral on the right hand side: 



JO 



\X ( t - t) h( U + U, V0) (t) || L cc cLt < c min{ 1 , t~ * }pS* ( t) 

+ cmin{l,rf}sJ x (OsJ(f), 



thanks to Lemma 5.2 of [2] . 

Now we need to estimate the contributions of the diffusive part of the hyperbolic Green function 

f \\K{t-r)h{JJ+U,V(l>)[j)\\L^dr 
Jo 

rt n 



J"f n rt n 

£ \\K lj + i{t -T)p{T)d Xj <p{T)\\ L °°dT + \ £ ||^JCi ; - + i(f-T)p«/>(T)||L<x.dT 
j=i JO j=i 

<Sj(r)s|^(?) f min{l,(f-T)-?"2}min{l,T"t}dT 

J 

+ cf min{l,(r-T)"? _1 }pe" &T rfT + c f minjl, (f-T) - 2 -1 }p||poll£idT 
Jo Jo 

< minjl, r*}sj(flsj x (fl + c(min{l, ri _1 }p + r*p||poll L i). 
In conclusion, we obtain 

II UCflOi- <c [min{l, r* } II t/ ll L i + e" ci BObllfl» + minjl, r^}E 2 s 
+ min{l, fq^(t)£ s+ i + min{l, rhpS^it) 
+ min{l, rf }sJ x CflsJ(f) +p(min{l, r f" 1 }p+ r?p||p || L i)] . 



n n 



Let us recall that S^, < C. Then, substituting inequalities (44) and (45) in the previous one and multi- 

n 

plying by max{l, 1 1 }, we obtain 

i?!(0<C(l + i?J(f)), 

where C = C(.E s+1 ,D s+1 ,p). 

In conclusion, if initial data and the constant state p are sufficiently small, for the L°°-norm of the solution 
{U,(p) we obtain the following estimates 

||Lf(f)|| L co< minjl, r^}C(E s+1 ,D s+1 ,p), 



||0tt)|| LOO < minjl, t hc{E s+1 ,D s+l ,p). 

□ 
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